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1. Introduction. We consider finite element approximations to the problem 

u t - Au = / in fi x (0,T], 

u = 0on9ftx[0,T], (1.1) 
u(x, 0) = Uq(x). 

Here ft C R™ (n = 2, 3) is a bounded polyhedral domain, / is sufficiently smooth, Uq 6 
Loo(^), and u is a weak solution to ([IT]) lying in L oo (0, T; L oo (fi))niJ 1 (0, T; 
(a subset of C*°(0,T; L^). 

Adaptive finite element methods for approximating solutions to parabolic partial 
differential equations are popular because of their ability to efficiently resolve singular- 
ities and other rapid local variations in solutions. While most adaptive finite element 
methods are designed to control only energy norms of solutions, in many applied 
problems the goal output of a finite element computation is related to some other 
norm or functional of the solution. In this work we address control of the maximum 
error \\u — w/JI^^x^t]) f° r nn he element approximations of u. Ensuring good 
pointwise approximation of u is natural in many situations where u represents some 
physical quantity. Pointwise error control is also a natural goal when computing free 
boundaries, for example via level set methods (cf. [DDE05 ). Several recent papers 
have addressed adaptive finite element methods for controlling pointwise errors in el- 
liptic problems (cf. |Noc95j . [DDPOOj . [NSV03j . |NSV05j . |NSSV06j . [De06] , [De07] h 
However, the only previous pointwise a posteriori estimates for parabolic problems 
that we are aware of are contained in |EJ95j and [BoOOj , which we describe below. 

The goal of this work is to prove a posteriori error estimates in £oo(£! x [0, T]) for 
semi- and fully-discrete finite element approximations to (|1.1[) . For practical purposes, 
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our main result is an easily-computable error estimator for the backward Euler finite 
element discretization of (jl.ip . In order to describe this estimate, we introduce some 
definitions and notation. Let = to < t\ < ... < = T, I{ = (£j_i,£j), and 
Ti = ti — U-%. For each < i < N, let % be a triangulation of fi. We place only 
standard restrictions on the triangulations, requiring in particular that all triangles 
have aspect ratios that are uniformly bounded with respect to i = 0, . . . , N and that 
the triangulations are "edge-to-edge" (i.e., hanging nodes are not allowed). Let Sq 
be a finite element space consisting of the functions that are continuous piecewise 
polynomials of degree k on % and which are on dil. Letting v l (x) = v(ti,x) for 
any function v defined on ft X [0,T], we discretize the weak form of (11.11) by letting 
m° £ Sq approximate u and then defining u l h 6 Sq, 1 < i < N, via the implicit Euler 
recursion 

- / (u£ - ui -1 )& da: + f V<V0idx= f f&dx for all cf>i £ Sq. (1.2) 

In addition, let g l = f l — Uh " h — , i > 1. The definition of g° is slightly different; cf. 
§4.2. Then for any 1 < j < N, 

\\ uJ - Khootp) <ho - WhIUcCn) 

+ C (n) (2 + c(n) In (In h) 2 max f^oK,^) 

Tj 0<i<j ( 13 ) 

+ E / ll/-/*IU.(n)dt+5l|ff < -fl < - 1 |U oo(n ). 

i=l 

Here c(n) = 2 „^ 2+1 , ^00,0(^15') is a standard and easily-computable residual-type 
estimator depending only on v,h, g l , and the mesh 7^ , and h is the minimum diameter 
of elements lying in U^L^. A more precise definition of £00,0 is provided in W2.31 and 
Theorem 14.71 contains a precise statement of results. 

The error bound on the right hand side of (|1.3|) consists of: 

• an initial data estimator \\uo — ^lli^n)! 

• a spatial estimator C(f2)(2 + c(n) In ^-)(lnh) 2 maxo<i<j £00. o( u h, 9 % ) which ac- 
counts for spatial errors; and 

• a time estimator £Li Si, 11/ ~ fh^Q) dt + f \\g l - .g'" 1 || Loo (o) which ac- 
counts temporal errors. 

Note that the constants in front of the initial data and time estimators are 1, so 
that estimating these terms does not require estimating unknown constants. While 
the spatial estimator does contain the constant C(f2), we shall show later that (7(f2) 
depends only on properties of an underlying elliptic a posteriori error estimator. 

In addition to the bound (|1.3p . we also establish several other a posteriori es- 
timates for semidiscrete finite element approximations of u as well as for the fully 
discrete scheme (|1.2|) . While more difficult than (|1.3p to employ in practical adaptive 
codes, these estimates provide additional insight into a posteriori theory for pointwise 
norms, for example by establishing that the spatial terms in the middle line of (|1.3p 
can sometimes be bounded instead in a weaker negative norm. 

We next briefly compare the estimate ()1.3j) with the results of the previously cited 
works [EJ95] and |Bo00j . [EJ95] contains pointwise a posteriori estimates for parabolic 
problems discretized in space using standard piecewise linear finite element schemes 
and in time using a discontinuous Galerkin approach. However, these estimates are 
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stated without proof, do not appear to be proven elsewhere in the literature, and ad- 
ditionally assume restrictive hypotheses on the spatial finite element meshes. [BoOO 
similarly employs a discontinuous Galerkin time discretization and piecewise linear fi- 
nite element spatial discretization. The proofs in BoOO] involve proving quasi-optimal 
regularity estimates for a regularization of the parabolic Green's function, which is 
fairly difficult and also leads to an uncomputable a priori term in the upper bound. 
This method of proof essentially involves imitating in a parabolic context the maxi- 
mum norm estimates for elliptic problems originally proven in |Noc95] and DDPOOJ. 
We finally note that the results of both [BoOO] and [EJ95] are restricted to convex 
polyhedral domains. 

We emphasize some features of (|1.3[) that contrast positively with existing results. 
First, in the present work Q may be a nonconvex polyhedral domain (including a 
domain with cracks). In addition, we allow arbitrary orders of finite element spaces, 
and in particular do not restrict ourselves to piecewise linear elements. The estimate 
p.3p also does not require any impractical restrictions on the spatial mesh (in contrast 
to [EJ95] ). and does not contain any uncomputable terms depending on u in the upper 
bound (in contrast to [BoOOj. Finally, as we discuss further below, the proof of (|1.3[) 
is quite straightforward because we are able to reuse difficult elliptic results instead 
of imitating their proofs in a parabolic context. Thus the results that we present here 
are to our knowledge the first rigorously proven, fully a posteriori pointwise error 
estimates for finite element methods for parabolic problems. 

Essential to our development is the elliptic reconstruction technique introduced 
in [MN03] in the context of semidiscrete problems and extended to fully discrete 
problems in [LM06]. In essence, the elliptic construction TZut is a continuous elliptic 
representation of the discrete solution Uh, and Uh is the elliptic finite element approx- 
imation to IZuh with respect to the finite element space under consideration. Thus 
any a posteriori error estimates which are valid for elliptic problems on f2 may be 
used to estimate IZuh — Uh- The overall error u — Uh may then be bounded by first 
estimating u — Huh using PDE techniques for continuous parabolic problems and then 
estimating TZuh — u% using elliptic a posteriori estimators. The elliptic reconstruction 
may thus be regarded as an a posteriori analogue to the Ritz- Wheeler projection in 
standard a priori error analysis for parabolic problems (cf. |Wh73| . [Tn97j). 

Our use of the elliptic reconstruction technique in the context of pointwise error 
estimation for parabolic problems highlights its ability to fully leverage existing elliptic 
estimates. In particular, establishing a rigorous theory for a posteriori estimation 
of pointwise errors for Poisson's problem on polyhedral domains was a technically 
difficult enterprise (cf |Noc95j . [DDPOOj . [NSSV06] ). Relying on these elliptic results 
instead of mimicking them, our proofs employ only basic estimates for the heat kernel 
for the continuous problem and are relatively short and straightforward. A direct 
"parabolic" approach to the problem which does not use the elliptic reconstruction is 
also possible (cf [BoOO. ) , but such an approach is much more technically involved and 
as already mentioned has not led to the sharpest possible results. 

An outline of the paper follows. In §2, we provide common preliminaries and recall 
some facts concerning residual- type a posteriori error estimation for elliptic problems. 
In §3 we prove a posteriori estimates for semidiscrete approximations of while 
in §4 we consider a backward Euler time discretization of (jl.ip . 

2. Preliminaries. In this section we provide a number of preliminaries regarding 
heat kernel estimates, maximum norm a posteriori estimates for elliptic problems, 
and issues concerning mesh compatibility that arise in some of our estimates for fully 
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discrete schemes. 

2.1. Notation. We begin by defining suitable notation. W^(iY), where j is 
a nonnegative integer, will denote the standard Sobolev space of functions having j 

o 

derivatives in L p (fi), and W p (Q) will denote the functions in W p (Q) which in addition 
have zero trace on the boundary <9fL In addition, (•, •) denotes either the L2^l)- or 
[I/2(^)]"-inner product. Finally, we denote by L p ([a,b],W^(fl)) the functions whose 
spatial W-j norm lies in L p over the time interval [a, 6]. 

2.2. Analytical preliminaries. Our analysis relies heavily on properties of the 
heat kernel for the problem We sum up the necessary results in the following 
lemma. 

Lemma 2.1. Let Q C K" be a bounded open domain. Then there exists a Green's 
function G(x,t;y, s) for the problem That is, there exists a kernel G such that 

for u satisfying and (x,t) G SI X (0,T], 

u{x,t)= f G(x,t;y,0)u a (y)d y + f f G(x,t;y,s)f(y,s) dx dt (2.1) 
Jn Jo Jn 

is a weak solution of For s < t, G additionally satisfies the bound 

\\G(x,t;;s)\\ Ll(Q )<l. (2.2) 

Let also 2 < p, q < oo with 

n 1 1 

2p + q < 2' 

o 

Then we have in addition that G(x,t;-,-) G L q > ([0, T], W p , (Cl)), where p' and q' are 
the conjugate exponents to p and q. Also, 

\\ G ( x ,~t;->-)\\L q ,([0,T],Wi,(ty) < Cp:?( T ): ( 2 - 3 ) 

where C p ^ q depends on p, q, |fi|, and T. In addition, G(x, t; •, s) G Hq(Q) for < s < 
t. 

Finally, there exists a constant c(n) depending only on the space dimension n 
such that for s < t, 

\\G s (x,f,;s)\\ Lim <^- s . (2.4) 

Here we use the notation G s {x,t; •, s) = J^G(x, t; •, s). 

Proof: The existence of a Green's function satisfying (|2.1[) and (|2.3p may be 
found in Theorem 6 (p. 657) and Theorem 9 (p. 671) of the fundamental work |Ar68j 
of Aronson. To prove l|2.2p . we note that 

< G(x, t;-,s)< T(x, t] y, s) = (4tt(* - s))-' l/2 e'^T . (2.5) 

That is, the heat kernel on is bounded pointwise by the heat kernel on R". This 
fact may be proven e.g. by combining Lemma 7 (p. 677) of |Ar68j with Corollary 
8.2 and Theorem 8.3 of TJaOO]. Inequality (|2.2[) then follows from the fact that 
J K n T(x, t; y, s) dy = 1 for t < s. 
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In order to prove (|2.4p . we apply Corollary 5 of |Dav97j with 5 = | and e = | to 
find that 

G s (t, x- y, s) < 2 n l 2 -\t - S )- 1 (4 7 r(i - , S ))-"/ 2 e-fe^. (2.6) 
Performing a change of variables and integrating over K" yields (12. 4p with c(n) = 



3" 



2 " /2+1 ' __ ^ 

Remark 2.2. One may take advantage of the bound (|2.5|) and (|2.6p in order to 

explicitly incorporate dissipation of the heat kernel into (|2.2p and (|2.6p . For example, 

for n = 2 and s < t one may compute that 

||G(a:,t;-,s)|| Ll( n) < l-e'^r, (2.7) 
l|G s (x,i:,-, S )|U l(0 ) <^(l-e-^#). (2.8) 



It is possible to incorporate (12. 7p and (|2.8|) into pointwise a posteriori bounds with 
no great difficulty, and we shall pursue this briefly in Corollary 14.41 

2.3. Elliptic a posteriori estimates. In this section we cite several results 
that will enable us to bound a posteriori the elliptic reconstruction terms appearing 
in our estimates for parabolic problems. In this subsection we assume that v satisfies 

—At; = g in CI, 
v = on dCl 

where CI C K.", n = 2, 3 is a polyhedral domain. We additionally assume that T is a 
shape-regular simplicial decomposition of Cl, and define the Lagrange finite element 
space S = {wh G ff 1 (il) : wu\k is a polynomial of degree k on K, K 6 T}. Let also 
Sq = S DHq^CI). Finally, let Vh £ Sq be the finite element approximation to v defined 

by 

VvhVwh dx = fwh dx for all Wh 6 So- 



Our parabolic results assume a posteriori bounds for \\v — v^Wl (f2)j an d in some 

i + i 

p p 



circumstances also for \\v — Vh\\ w -^^y With 1 + ^ = 1, here 



IMIw- 1 ^) =sup{(w,x) : z S Wl,(n),\\z\\ w i /{n) = 1}. 

We shall employ residual-type estimates. We first define the jump residual [Vv/j] on 
an (n — l)-dimensional element face e = K\ n K2, where K\, K2 £ T. Let ft be an 
arbitrary unit normal vector on e, and for x G e let 

[Vwfe](a;) = lim(Vu/,.(a; + 5ft) — Vvh(x — 5n)) ■ n. 

8— >0 

Let also be the diameter of the element K. For 1 < p < 00 and j > 0, we then 
define the elcmcntwise error indicator 

V P ,-j(K) = h 2 +i\\g + Av h \\ Lp{K) +h£ 1+lp mv h }\\ Lp{9K) . 
Finally, we define the global estimator 

^, 9 )4 (E ^*- ( ^r i£p<00 - (2.9) 

I max^grTyoo -j{K), p = oo. 
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We first quote an elliptic pointwise estimate which holds for all orders of finite 
element spaces and all polyhedral domains; cf. [NSSV06] for a proof. 

Lemma 2.3. Assume that f2 is an arbitrary polyhedral domain in W l , n = 2, 3, 
and let h = mixiK^T hn . Then 

When allowed by the domain geometry, it may be advantageous to instead mea- 
sure the error in a negative norm. In such cases we shall employ the following result. 

Lemma 2.4. Assume that Q is a convex polyhedral domain in W 1 , n = 2,3, where 
the maximum vertex opening angle (for n = 2) or edge opening angle (for n = 3) is 
denoted by to — —, a > 1. Assume also that the degree of the finite element space Sq 
is at least 2, that is, So contains the continuous piecewise quadratic functions. Then 

■f or dh" < P < °°> 

\\v-vh\\ w -i m < C(p,Q)£ p ,-i (2.10) 

Proof: Assume that —Aw = z, where z 6 W^i(Cl) for | + ^ = 1. Combining the 
comments of §4.c and Corollary 3.9 of |Da92] yields the regularity result 

IIHIw^n) - c (p'^)\\ z \\w^(u)- 

Given this estimate, the result (|2.10p may be obtained using a duality argument and 
standard techniques for proving residual-type a posteriori bounds. □ 
Remark 2. 5. We emphasize that Lemma 12.41 only holds on convex polyhedral 
domains. It should be possible to similarly prove usable negative norm estimates on 
nonconvex polygonal domains in R 2 , but explicit information about corner singular- 
ities would appear in such estimators (cf. [LN03 for analogous global L 2 bounds on 
nonconvex polygonal domains) . On nonconvex polyhedral domains in M 3 , such a re- 
sult is much less practical as the precise nature of vertex singularities is often difficult 
to ascertain. 

Remark 2.6. Care must be taken when employing (|2.10p as the constant C(j>, ft) 
degenerates as p — » oo, and possibly as p approaches the lower bound —35- as well. In 
particular, one is not able to choose p(a) so that C{p{a), Q) remains bounded as the 
maximum edge opening angle approaches 7r (i.e., as a — > 1). 

2.4. Compatible meshes and estimates for differences in finite element 
solutions. Some of our fully discrete a posteriori estimates require bounding elliptic 
finite element errors of the form v± — 1>2 — {vhi — ^2), where Vhi and Vh2 he in different 
finite element spaces (in particular, in finite element spaces defined on meshes at 
adjacent time steps). In this subsection we make assumptions and definitions on the 
pair of meshes that allow us to establish such estimates. We follow closely Appendices 
A and B of |LM06j . so we shall be brief and refer the reader to that work for more 
details. 

Two simplicial decompositions T± and T2 of Q are said to be compatible if they are 
derived from the same macro triangulation M by an admissible refinement procedure 
which preserves shape regularity and assures that for any elements K G T\ and K' £ 
T2, either KDK' = 0, K C K', or K' C K. The bisection-based refinement procedure 
used for example in the ALBERTA finite element toolbox (cf. [SS05]) is known to be 
admissible. 
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There is a natural partial ordering of compatible triangulations, with T\ < T 2 if 
I2 is a refinement of T\ . The finest common coarsening Ti A T 2 of 7[ and T 2 is defined 
in a natural way, and h — max(/ii, h 2 ), where hi, h 2 , and h are the local mesh size 
functions for 7i , T 2 , and 7{ A T 2 . Finally, let Si and S 2 be finite element spaces of 
degree k on 71 and 72. 5 = Si (1 S2 is then the corresponding space of degree k on 
71 A72. Standard interpolation inequalities hold for all of the above-mentioned spaces, 
though the constants in these bounds may depend on the number of refinement steps 
used to pass from 71 to T 2 . We again refer to Appendix B of [LM06 , Appendix B for 
more discussion. 

Let Vhi £ Si and Vh2 G S2 be the elliptic finite element approximations to v\ 
and v 2 lying in finite element spaces Si and S2 defined on different meshes T\ and 72 , 
respectively. Here we assume that — A«i = gi, — £±v 2 = g 2l and v± = v 2 = on d£l. In 
essence, Lemma 2.3 and Lemma 2.4 still hold, but with the local mesh size h replaced 
by the local mesh size h of the finest common coarsening. Let S i; i = 1,2, be the 
union of the faces of elements lying in %. For K E Ti f\T 2 , let = (Si U S 2 ) n K, 
where K is taken to be closed. 

For 1 < p < 00 and j ; > 0, we then define the elementwise error indicator 

= h 2 + j \\ 9l - .92 + Av hl - Av h2 \\ Lp{k) + ^J 1+ ^||[VKi - « M )]|| ip(Ejt ). 
Finally, we define the global estimator 

P (v v o a r r ^-/ ( ^eT 1 AT 2 ^(W 1/ M<p<^ 

c- Pt -j(v h i - v h2 ,gi— g 2 , J 1 , J 2 ) — < „ /r >» 

[max^^^oo-jtiT), p = oo. 

Corollary 2.7. Assume that fi is an arbitrary polyhedral domain in M™, n = 
2, 3, and that 7{ and T 2 are compatible triangulations. In addition, define h = 
min^ga min(ft,i(x), h 2 {x)). Then 

hi - V2 - (vhi - Wft2)|U (X> (n) < C(^)(lnA) 2 ^oo,o(w/ii - v h2 ,g! - g 2 ;T 1} T 2 ). (2.11) 

Here C(f2) depends on the number of refinement steps used to pass from T\ to T 2 . 

COROLLARY 2.8. Assume that Q is a convex polyhedral domain in R ra , n = 2, 3 ; 
where the maximum vertex opening angle (for n = 2) or edge opening angle (for n = 3 ) 
is denoted by u> = — , a > 1. Assume also that the degree k of the finite element 
spaces Si and S 2 is at least 2, that is, both spaces contain the continuous piecewise 
quadratic functions. Assume also that T\ and T 2 are compatible triangulations. Then 

f° r S^i < P < 00 ' 

\\vi -V2- (vhi - Vh2)\\ w -i(n) - C(p, fi)£p,-i(vhi - Vh2,gi - 52; 7i,T 2 ). (2.12) 

Here C(p, £1) depends on p and the number of refinement steps used to pass from T\ 
to T 2 . 

The main observation used to derive the above corollaries is the fact that (v\ — 
v 2 ) — (vhi — Vh2) is Galerkin orthogonal to the space Si n S2. The compatibility of T\ 
and 72 ensures that this intersection is rich enough to obtain (|2.11|) and (|2.12|) . The 
proofs otherwise follow closely those of Lemma 12.31 and Lemma 12.41 and we do not 
give details here. 

3. Analysis of the semidiscrete scheme. 
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3.1. Semidiscrete finite element approximation. For simplicity in handling 
finite element approximations, we assume that f2 is a polyhedral domain in R™, n = 
2, 3, and that 7^ is a simplicial decomposition of fl. We emphasize that we admit 
here general polyhedral domains, including those having cuts or cracks (as described 
in the standard reference |Da88j ). Let S C if 1 (SI) and So = S D Hq(H) be standard 
simplicial Lagrange finite element spaces as in §2.31 The semidiscrete approximation 
Uh G C([0, T], So) of u then satisfies 

(uh,t,v h ) + {Vu h , Vv h ) =(/, v h ), v h G S and t G (0, T], 

Mo) =< 

where G So is a finite element approximation to uq. We let Ph : L>2 — > S be the 
Li projection onto the finite element space S, and additionally define the modified 
discrete Laplacian -A h (t) : H^(Q) -> (S + Phf{t)) by 

(Vtt, Vv h ) = (-A h (t)u, v h ), v h G S . (3.2) 

From ()3.1|) and (|3.2j) . we have the pointwise formula 

- A h u h = P h f - u h>t . (3.3) 

Remark 3.1. Our definition of —Ah is nonstandard in that here — AhU has nonzero 
boundary values which depend on the data / in addition to u. We use this definition 
in order to maintain consistency in the pointwise relationship (13. 3|) . In particular, 
because Uh,t = on <9f2, we also have Phf = —Auh on Oil. Note that we instead 
could enforce this relationship by letting — AhUh G So and taking the L2 projection 
of / onto So- This distinction will make little practical difference in our development, 
but it is possible to define the elliptic reconstruction in such a way that — AhUh and 
Phf must be computed a posteriori (cf. |LM06] L 

3.2. Elliptic reconstruction for the semidiscrete problem. Given a finite 
element approximation Uh, we define its elliptic reconstruction TZuh G Hq(H) by 

(VKuh, Vv) = (g,v),ve ff x (O), (3.4) 

where g = / - u h .t- We thus have -A h u h = Phf - Uh,t, -AlZu h = / - u h .t, and 
—Au — f — u t . Note that TZ = lZ(t) is a time-dependent operator, but we shall depress 
its dependence on t in the sequel in order to avoid unnecessary clutter in our notation. 
We will deal more explicitly with the time dependence of the reconstruction operator 
in our analysis of the fully discrete scheme. Also note that we may differentiate (|3.4[) 
with respect to t to obtain 

{V{nu h ) t ,Vv) = {g u v), 

where g t = ft - u h<tt . 

Combining (|3.2p . (|3.3p . and (|3.4p . we find that Uh and IZuh satisfy the Galerkin 
orthogonality relationship 

{V{TZu h - u h ), Vv h ) = 0, v h G S . (3.5) 

Differentiating (|3.5[) with respect to t also yields the time-differentiated Galerkin or- 
thogonality relationship 



(X7(Ku h - u h )t, Vv h ) = 0, Vh G So(fi). 
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In addition, it is easy to calculate that for < t < T and v £ H ( 



(n), 



((it - TZu h ) t ,v) + (V(u - Tluh), Vu) = ((uh - fcu h )uv). 



(3.6) 



The fact that u — IZiih thus satisfies a parabolic equation will play a fundamental role 
in our development. 

Remark 3.2. In MN03], the elliptic reconstruction is defined by 



From (|3.3|) , we see that —A^Uh + f — P^f = f — Uh,t so that (|3.T|) and (|3.4[) are 
equivalent. In fact, the elliptic reconstruction allows us to write the pointwise form 
(|3.3p of the discrete equation as 



The above equation does not involve the discrete Laplacian and thus allows for a 
straightforward comparison with the PDE leading in its weak form to (|3.6p . 

We use the definition (|3.4[) because in what follows we employ residual estimators to 
estimate the elliptic error TZuh — Uh- These estimators require pointwise access to the 
right-hand-side data for TZuu- It is not practical to directly compute —AhUh or Phf, 
and in MN03] , the authors develop an expression for the residual that does not involve 
the operators — or Ph (cf. p. 1592). Thus the definition in |MN03 emphasizes the 
underlying structure of the reconstruction operator, that is, 1Z = (— A) _1 (— A/j) up 
to terms that are /^-orthogonal to the finite element space. Our equivalent definition 
instead reflects the practical concern of computing using the resulting a posteriori 
error estimates. 

3.3. Semidiscrete reconstruction results. Our a posteriori estimates are 
based upon the following theorem. 

Theorem 3.3. Let the assumptions and definitions of §3.1 and §3.2 hold. Then 
for any < to < T , 



(VKu h ,Vv) = (-A h u h + f- P h f,v), v e H*(Q). 



(3.7) 



Uh,t - ATZu h = f. 




(3.8) 



\{u-u h )(t )\\ Lx (n) <\\{R>u h - u h ) (to) Hz^n) + \\{u-Ku h )(0)\\ Lx{n) 
+ C p , q {to)\\{u h - 'R-u h ) t \\ Lq{[a to] . w -i {n)) . 



(3.9) 



Proof: For any Xq €: fl 



\(u - Uh)(x ,t )\ < \(u ~ Ku h )(x ,t )\ + \\(TZu h - u h )(to)\\ Loo (n). 



Using dUTJ and ([3T6]) . we find that 




+ 




G(x ,t ; y, s)(u h - Ku h ) t (y, s) ds 
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The first term on the right hand side above may be bounded by combining Holder's 
inequality with (|2.2[) . In order to bound the second term, we apply (|2.2p to obtain 
(pTSj) or (E3D to obtain □ 
Remark 3.4- We have assumed a polyhedral domain and a specific type of finite 
element space in Theorem 13.31 but similar results hold under more general circum- 
stances. Indeed, the analytical results of Lemma l2~Tl hold on general bounded domains, 
and only those estimates along with the relationship (|3.6|) are used in the proof. The 
bound (|2.2p also holds for a fairly general class of elliptic differential operators, though 
in more general cases one must perhaps replace 1 on the right hand side by an un- 
known constant with unknown dependence on T. The reconstruction technique is 
thus able to transfer most difficulties and issues concerning the precise nature of the 
finite element approximation (including for example the type of elements used and 
difficulties arising from finite element approximations on nonpolygonal domains) to 
the a posteriori estimation of elliptic errors. 

3.4. A posteriori error estimates for the semidiscrete problem. In this 
subsection we shall estimate the right hand sides of (|3.8[) and (|3.9p a posteriori using 
the residual estimators of §2.31 We first present an estimate which is valid for general 
polyhedral domains. 

Theorem 3.5. Let Q C K™, n = 2, 3, be an arbitrary polyhedral domain, and let 
Uh G So be a standard Lagrange finite element approximation defined on an arbitrary 
shape-regular simplicial decomposition of having minimum mesh diameter h. Then 
for < t < T, 

||0 - u h )(t )\\ Lgo{n) <\\u Q - ul\\ Loo( Q) 

+ C(Cl) (In h) 2 [8^ , (u h (0) , g(0) ) + £, x .o (u h (t ) , g(t ) ) 

+ l|£oo,o(wM'fft)ll-Li((0,T))] ; 

where -Eqo.o * s the L^-type residual estimator defined in 12.9\) . 

Proof: We proceed by bounding the terms in (|3.9p using Lemma T2.3I Recalling 
(pTl)) and (f3~5]) . we find that 

\\(TZu h - u h ){t Q )\\ Lao{n) <C(n)(ln£) 2 £ 00 ,o(ufc(t),5(t)) 

and similarly 

\\nu h (0) - u°\\ Lx{n) <\\u - u° h \\ Loom + C(n)(lnfe) 2 £ oo>0 «,fl(0)). 

Finally, 

\\(u h - TZu h ) t \\ Ll{[0it]MQ)) < C(n)(lnft) 2 [|5 00)0 (it/ l ,t,fl'i)||i 1 ([o ) t]). 

Inserting the above inequalities into p.8p completes the proof. □ 
We next present a theorem which allows us to bound the main error term with a 
higher-order estimator. However, this estimate only holds for quadratic and higher- 
order elements and convex polygonal domains. 

Theorem 3.6. Assume that Q is a convex polyhedral domain in W l , n = 2,3, 
with maximum vertex (if n = 2) or edge (if n = 3) opening angle — , a > 1. Let also 
2 < P, q < oo satisfy ^ + ^ < | and < p. Finally, assume that the polynomial 
degree of the finite element space So is at least two. Then with h = min/f e r h diam(X), 

\\(u-u h ){to)\\ Loo (n) < \\u -u° h \\ Loom 

+ C{n)(lnh) 2 [£ oofi (u h (0),g{0))+£ ODfi (u h (to),g{to))} 

+ Cp, g (t )C(p, ti)\\£p,-l(Uh,t,9t)\\L g ([0,to])- 
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Proof: The proof is completely analogous to that of Theorem 13.51 above, the only 
difference being that we now employ Lemma HOI in addition to Lemma |2~3"1 □ 



4. Analysis of the fully discrete scheme. 

4.1. Fully discrete finite element approximation. As in EJTJ let = to < 

t\ < ... < t N = T, Ii = (£,•_!,£,•), and 73 = |/,|. For each < i < N, let % be a shape- 
regular simplicial decomposition of fl. Let also S l be a space of continuous piecewise 
polynomials of degree k on %, k > 1, and Sq — S l (~l Hq(SY). Defining v l (x) — v{U, x) 
for v £ C(fi x [0, T]), we discretize the weak form of (jl.ip using the backward Eulcr 
method as follows. Let w° £ Sq be some approximation to u l h £ Sq, 1 < i < N, 
is then defined via the recursion 

( u i - u^ 1 \ 

-* & + (V<, V&) = (/\ for all 0, G S*. (4.1) 



In order to obtain a discrete approximation to u on the whole parabolic domain 
Q x [0,T], we interpolate the functions u % h linearly between and tf. 

u h { X .t) = (1- ~ ^v rV) + < * < u. 

n n 

Finally, we define u\ t = Jjtiftli^, that is, 

<M = <^~<^\ i > 1. (4.2) 

n 

Next we define L2 projections onto S l and Sq. For < i < N, we define P^ : 
L 2 {Q) -> S 4 and i*„ : L 2 (fi) S^ by 

(P^u,Vi) = (u,Vi) V u, € S l , 

The discrete Laplacian — Al : Pq (CI) — > Sq + P/j/ z is then given by 

(-A>, ^) = (Vu, V«i) for all v t £ Sj. (4.3) 

The weak-form fully discrete scheme (|4.ip may easily be transformed into the 
pointwise equation 



pi 



n 

Referring to (|4.2p . we thus find that 



*Wh = nr 



pi 7; 4 -! _ 



A?X = P^ 1 + i > 1. (4.4) 



4.2. Elliptic reconstruction. We define the elliptic reconstruction by first 
defining it at the time nodes and then interpolating linearly between them. Following 
(HH) , for < i < N we let 

(V7T<, V«) = (g» V v G tf^n), (4.5) 
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where 

We then obtain the time-continuous elliptic reconstruction 

Ku h = (1 - t —^)TZ l - 1 u\- 1 + t—^Wul, U-i <t<U. 
Using the formula (|4.4[) , it is easy to compute that 

r - 4 )t = [-am + ir - pin + i ul ~ PlfiUh ]• (4.7) 

n 

This relationship combined with the definitions (|4.3|) and (|4.6|) yields the Galerkin 
orthogonality relationship 

(V(< - KVJ, Vw h ) = V u h e Si (4.8) 
Similarly, we have on Ii that 

(v«, - (Hu h ) t ), Vu h ) = v u fe e s^ 1 n s**. (4.9) 

We finally state an error equation which will play a fundamental role in our 
analysis. For t 6 Ii and (f> £ -^o(^)> 



((u - 7fc*/,)t> 0) + (V(u - ftu A )> V0) =((« h - ftu h )t, 0) + (/ - / 



(4.10) 



Comparing (|4.10[) with (|3.6p . we see that (|4.10[) contains additional terms which all 
result from the time discretization of (jl.ip . 

Remark J^.l. From (|4. T[) . we see that the elliptic reconstruction for the fully 
discrete problem lifts the sum of a discrete Laplacian term, a spatial data approxima- 
tion term, and a mesh coarsening term. Note that the last term is nonzero only if 
So -1 Sq, that is, if the mesh is coarsened in the i-th time step. Our definitions of TV 
for i — and i > therefore differ only by the exclusion of the mesh coarsening term 
when i — 0. The definition of 1Z L in LM06J does not incorporate the data approxi- 
mation and mesh coarsening terms in (|4.7[) into the right hand side of the equation 
solved by l? 1 ^. The advantage of including these terms is that the error equation 
(|4.10[) now includes only the term ((ith — TZuhjti 4>) P ms terms which result from the 
discretization in time, resulting finally in a posteriori estimators which have a simpler 
structure. 

4.3. Fully discrete reconstruction results. Here we present three alternative 
results. 

Theorem 4.2. Let all assumptions and definitions be as in §^.J and $4-2. Then 
for any 1 < j < N , 

\\(u - u A )(*,-)IU=o<n) <\\(Ku h - UfcXtjOIUoopi) + ||(« - Ku h )(0)\\ Loo{n) 

+ \\(u h - 'R-Uh)t\\Li([0,t j ];L oo (n)) 

i . ( 4 - U ) 



i=l Jl i 
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Alternatively, for 2 < p, q < oo satisfying ^ + | < i and /or any i < j < N , we have 

\\(u - u h )(t 3 )\\ Loa(n) <\\(Ku h - u h ){t 3 )\\ Loa{n) + \\{u- Ku h ){Q)\\ Lxi{n) 
+ C P , q (tj)\\(u h - Ku h )t\\ Lq{[o t]] . w -i {n)) 



(4.12) 



Finally let c{n) be the constant from t\2.4l ■ Then 



(u-u/,)(tj)llioc(n) < Hwo-^olU^co) 

+ (2 + c(n) In ^) max. ||7T< - <|| Loo( n) 



E / ll/-/il^(n)dt+||| 5 



5 l II £«,(«)■ 



Proof: We first fix a point x G 51 with \\(u - = l( u - u h){xo,t 3 ) 

and compute 

(w - u/,)(;ro,ij)| < |(" - Uu h )(xo,tj)\ + \\(TZu h - Wft)(tj)IUoc(n)- 
Using (H3J), (jlTTUl) , and we find that 

(u-Kuhjixctj) = / G(x Q ,t J - 1 y,0)(u-TZu h )(y,0)dy 



3 J / u l*«' 1 )' 

Jn 

nG(x ,tj;y,s)(u h - TZu h ) s (y,s) dy ds 
.1 

E/ / G(xo,t J - 1 y,s)(f-f)(y,s)dyds 
i=1 Jii Jn 

E / / G(x Q ,t J - 1 y, S )(l--^^)(g*-g>- 1 )dyd S 



i=i 



(4.14) 



Jn 

<\\(u-Ku h )(0)\\ Loo{n) + / / G(x Q ,t j ;y,s)(u h -Ku h ) s {y,s)dyds 

Jo Jn 



E / ii/-/iu»co)d*+f lis -//' 1 Lol >, 

i=l 



The term JL J J„ G(a/o, tj-; y, s)(uh — TZu^tiy, s) dy ds above may then be bounded by 
using (|2.2p in order to obtain (|4.1ip . or by using (|2.3D to obtain (|4.12|) . 

In order to prove (|4. 13|) . we begin as in (|4. 14|) . split the second integral into two 
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integrals, and perform integration by parts in time on [0, tj-_i] to compute 

(u-TZuh)(x ,tj) = [ / G(x ,tj;y,0)(u - TZu h )(y,0) dy] 
Jn 

+ [/ / G(x ,tj;y,s)(uh-1lu h ) s (y,s)dyds] 
Jo Jn 

+ E/ f G(x ,t j ;y,s)(f-f)(y,s)dyds] (4-15) 
i=1 J/, Jn 

+ E/ / G(x ,^;y, S )(l-i^i)( 5 ' ; -.g i - 1 )dyd S ] 

=[/] + [/ j + jij + jv + v] + [vi] + [vii], 



where 



1=1 G(x ,tj;y,0)(u-TZu h )(y,Q)dy, 
J si 

11= / G(x ,tj;y,tj-i)(uh-lluh)(y,tj-i)dy, 



n 



III = - / G(x ,tj;y,0)(u h -1lu h ){y,0)dy, 

Jn 

IV = - I / G s (x ,tj;y,s)(u h - Tlu h ){y,s) dyds, 

Jo Jsi 

V= / G(xQ,t j ;y,s)(u h -Ku h ) s (y,s)dyds, 

Jtj-i Jn 

j 



VI = J2 / G(a;o ) t i ;i/,«)(/-/ i )(y,a)di/da ) 

i=l 

VII = Y,[ I G{x^t ]]y , S ){l- S -^-){^ -g'-^dyds. 
i=1 J/j T i 

In order to bound the terms in (|4.15|) . we use ()2.3|l to find that 



1 + 11 + 111= G(x o ,t J ;y,0)(u-u h )(y,0)dy 



+ / G(x ,tj;y,tj-i)(uh - Uuh){y,tj-x) dy 
Jn (4.16) 

<||G(a;o,*j; • ) 0)|| il( n)||wo - ttfcolU^n) 

+ l|G(xo,tj;-,ti-i)||L 1 (n)||(Mfe-^ h )(t i _i)|| Loo (Q) 

<||«o - u hQ \\ Loa{n) + \\{u h - 1lu h ){tj-i)[[ Loa{Q) . 
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Employing (|2.4[) . we obtain 

IV < \\G s (xo,t ;-,s)\\ Ll{n) \\(u h -Jiu h )(s)\\ Loa{n) ds 
Jo 

c ( n ) 

<\\uh — T^UhW Z.oo(f2x(0,i3-i)) / ds (4.17) 

Jo l j ~ s 

< c ( n ) l n ^ o< max 1 1|< - R i <u i h \\ LoB (p)- 
In order to bound the term V , we calculate 

V <||G(x ,to,-,-)llLi(Ox/ J )ll('"/ l - ^/ l )t||i 00 (f2x/ J ) 

< T ill " lUoo(O) 

<K - Wv? h \\ Loom + IK' 1 n^ui-'U^ny 
Finally, we compute directly that 



i=i ■'-f* z 



^Ik^n). (4.18) 



Collecting the previous inequalities and inserting them into (|4. 15|) completes the proof 
of gH3J) . □ 
Remark 4-3- One may approach the proofs of (|4.1ip and (|4.13p of Theorem 14.21 
and l|3.8p of Theorem 13.31 from the perspective of semigroup theory instead of using 
fundamental solutions. Let E be the semigroup generated by the Laplace operator, 
i.e., let E(t)vo be the solution of v t ~ Av — 0, v(0) — vq. Then using Duhamel's 
principle we have 

u(t) = E(t)u + [ E(t- s)f{s) ds. 
Jo 

Assume that the stability and strong stability properties 

\\E(t)\\L a> ^L BB <l,t>0, (4.19) 
\\E'(t)\\ L ^ Lao <j,t>0 (4.20) 

hold. The bounds (13. 8p and (|4.13p are then easily obtained by respectively combining 
(f3~6]) and ([4~T0]) with (|4~19|) . (ILLS]) may be obtained by combining (|4l"9|) . (fL20|) . and 
(I4.10p . We note, however, that we are not aware of a proof for the results (|4.19p and 
(|4.20p under the weak restrictions that we have placed on fl. The standard reference 
[St 74 , for example, assumes that dfl is C 2 in order to obtain the analyticity of E in 
C° and thus obtain P~2Tj]) . 

We finally note that it is possible to sharpen the reconstruction estimate (|4.13p 
somewhat, though at the cost of a more complex result. In particular, we may employ 
the estimates ()2.7p and (|2.8p that reflect the dissipation of the heat kernel when n = 2 
and also accumulate the spatial errors in t\ instead of l^. Let 4>i{s) be the piecewise 
linear "hat" function satisfying (j>i(ti) — 1, 4>i{t m ) = for m ^ i. Instead of (|4.17p . 
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we may then calculate 



IV < \\G s (xo,to;;s)\\L^Q)\\(u h -1lu h )(s)\\ Lao ^)^s 
Jo 

<c(n)[\\u° h -n u° h || ico(n) / Msh (1-e ^T^)ds 

Jo l j — s 

3 fti + l 1 diam(i"2) 2 

+ E 114-^411^(0) / M°h (1-e W^) da 

i=i Jti-i l o - 8 

+ \\u]- 1 -W- 1 u]- 1 \\ Loom / ^._ l(s ) (l- e ^5 ) ds ] (4.21) 



tj-2 Z 3 S 



/y-^ diam(f2) 



<cW[ 2fe _ fi) (i - e ^-'i' )||< - n u° h \\ Loo(n) 

7—2 o 

i=1 l i j z i i j z i+\ 

diam (f2 ) 



2t 3 



^ )ii«r 1 -^ i u r i iiL .(a ) ]. 



Employing (|4.2ip instead of (|4. 17)) and similarly inserting IpTF)) into (|4.16p and (|4. 18)) 
leads to the following result. 

Proposition 4.4. // the spatial dimension n = 2 and i/ie conditions of Theorem 
are met, then 

(u - u h )(tj)\\ Loo ^ < /5 4 (0)|[u - u h o\\ Loo (Q) 



Z Lj Zi lj 

+ (2 + cCnJ^-Afe-!))!!^- 1 - ^'-^ri^co) + 2||u> - WJ| ioo(0) 



*j - ti+i (4.22) 



i=l i=l 



Iff'-fl^lU.tn). 



Here 



di a m(n) 2 

/? 4 (s) = 1-e 4<t ^ 3) , 



/3 9 (s) = 1-e 9(t i" s) 



In contrast to (|4.13p , (|4.22p reflects damping of the error effects from times U 
for which tj — t{ >> 0. This damping is reflected in two ways. First, in (|4. 22)) the 
spatial indicators \\u l h — T^uJJI^^) accumulate in l\ instead of in 1^ as in (|4.13p . 
The weights i( t T \. + t .^j" 1 ) in (|4. 22)) deemphasize the spatial indicators for times 
ti « tj relative to contributions from times ti ss tj. This yields a sharper bound than 
(|4.13)) , where the spatial indicators are all weighted equally. Note that similar weights 
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can also be obtained when using a duality- reconstruction combination (cf. [LM07 ). In 
addition, the dissipation estimates (|2.7p and (|2.8p are used above to quantify damping 
in the initial data error (first line above) and time discretization indicators (last line) 
as well as in the spatial errors. While not dramatic, these dissipation effects can 
be substantial. If for example tj = diam(fi) = 1, then /^(O) « .222, corresponding 
roughly to one level of refinement in the initial mesh if piecewise linear elements are 
used. 

We finally note that (I4.22[) still holds when n > 2, but the dissipation weights 
Pi(s) and (3g(s) take a different form. 

4.4. A posteriori error estimates for the fully discrete problem. We 

finally obtain three different a posteriori estimates. 

Theorem 4.5. Let C K", n — 2,3, be an arbitrary polyhedral domain, and 
define h = mino<i<jv nunrgXi hx- Then for any 1 < j < N , 

\\(u - u h )(tj)\\ Loo (n) <IK - u° h \\ Loom 

+ C(n)(\nh) 2 [£ ooS) {v 1 h ,g j ) + 5 00 , «,<?o) 

3 

+ ^-.oK - U V . 9 l - <T *5 3-1, %)\ (4.23) 




Theorem 4.6. Assume that Q is a convex polyhedral domain in M. n , n = 2,3 
with maximum vertex (if n — 2) or edge (if n = 3) opening angle — , a > 1. Let also 
2 < P, q < oo satisfy + ^ < ^ cmd ^-j- < p. Finally, assume that the polynomial 
degree of the finite element space Sq is at least two. Then for 1 < j < iV, 

|[ («— tifc) (*,- ) | (n) < \\u - uIWl^q) 

+ C^Xln/iH^oK, P - u{ t ) + £^ Q (ulg Q )} 

( J \ V9 

+ C M (t j )C p (Sl) f^r-^^K-^- 1 ,^-^- 1 ;^,^) (4.24) 

+ E / iiz-ziiw^dt+^ii^-^iu^n). 



Theorem 4.7. Let Q C K", n = 2, 3, be an arbitrary polyhedral domain, and 
define h = mino<i<Ar nihiTgTi n T- Then for any 1 < j < N , 

\\(u - u h ){tj)\\ Loo{n) <\\u - ulh^n) 

+ (2 + c(n)ln^)C(ri)(\nh) 2 max ^(ul,^) 

Tj ' o<i<j ' ( 4 25 ) 

+ E / HZ-ZiUooWdt + lll^-^- 1 !!^^), 

1=1 jl * 

where c(n) — 2rl % +1 ■ 
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Proof: The proofs of the three preceding theorems follow easily by inserting the 
estimates of Lemma 12. 31 Lemma 12.41 Corollary 12. 71 and Corollary 12.81 into the esti- 
mates (|4~TTj) . P~T2")) , and P~T3]l while recalling the definitions ((43)1 and (|4~5)) and the 
Galerkin orthogonality results (|4.8|) and (|4.9p . □ 

Remark 4-8- In contrast to (|4.23|) and (|4.24[) . the estimate (|4. 25[) does not re- 
quire the computation of residual-based estimators with respect to a finest common 
coarsening. The fact that spatial errors in (|4.25p accumulate in L M in time as well 
as in space also is practically advantageous as it is easier to ensure that these errors 
are of the correct size at each time step in an adaptive algorithm. In addition, (|4.25p 
more readily lends itself to establishing a convenient bound for \\u — tt/iH^^xIQ,*-])- 
Thus for practical purposes, (|4.25|) is of greatest interest among the results in this 
paper. We also emphasize that (|4. 25(1 does not include any unknown time-dependent 
constants. In fact, all unknown constants in (|4.25| stem from the use of a posteriori 
error estimators for elliptic problems. 

While not as practically advantageous as Theorem l4.7i Theorem 14.51 and Theorem 
14 .61 also have interesting theoretical features. (|4.23|1 is interesting in that it includes no 
time-dependent constants of any sort. (|4.24ll bounds the spatial errors at intermediate 
times in a weaker negative norm in which u — Uh is of higher order for quadratic and 
higher-order finite element spaces. 

Remark 4-9. One may sharpen Theorem [477] as in Corollarv l4.4l The latter result 
may easily be adapted to obtain an a posteriori estimate, but the resulting estimator 
is cumbersome and we thus do not record it here. 
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